LOCAL MULTIGRID IN H(curl) 
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Abstract. We consider fJ"{curl, n)-elliptic variational problems on bounded Lipschitz polyhedra 
and their finite element Galerkin discretization by means of lowest order edge elements. We assume 
that the underlying tetrahedral mesh has been created by successive local mesh refinement, either 
by local uniform refinement with hanging nodes or bisection refinement. In this setting we develop 
a convergence theory for the the so-called local multigrid correction scheme with hybrid smoothing. 
We establish that its convergence rate is uniform with respect to the number of refinement steps. The 
proof relies on corresponding results for local multigrid in a (Q)-context along with local discrete 
Helmholtz-type decompositions of the edge element space. 
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1. Introduction. On a polyhedron fiC M'^, scaled such that diam(f7) = 1, we 
consider the variational problem: seek u G Hpo (curl, ft) such that 

(1.1) (curl u, curl v)^2(jj) + (u,v)^2(o) (f:V)^2(j^) Vv e ffpo (curl, f]) . 
" V ' 

= :a(u,v) 

For the Hilbert space of square integrable vector fields with square integrable curl 
and vanishing tangential components on To we use the symbol -Ffro (curl, $7), see [24, 
Ch. 1] for details. The source term f in (|l.ip is a vector field in (L^(17))^. The left hand 
side of (jl.ip agrees with the inner product of Hr^ (curl, f2) and will be abbreviated 
by a(u,v) ("energy inner product"). 

Further, To denotes the part of the boundary dfl on which homogeneous Dirichlct 
boundary conditions in the form of vanishing tangential traces of u are imposed. The 
geometry of the Dirichlet boundary part F^i is supposed to be simple in the following 
sense: for each connected component F^ of F^j we can find an open Lipschitz domain 
D.i C such that 

(1.2) Tiir\Ti = ri, a, nr2 = 0, 

and Qi and Qj have positive distance for i ^ j. Further, the interior of £7 U fii U 
^2 ■ ■ ■ is expected to be a Lipschitz-domain, too (see Fig. 15.21) . This is not a severe 
restriction, because variational problems related to (jl.ip usually arise in quasi-static 
electromagnetic modelling, where simple geometries are common. Of course, Fd = 
is admitted. 

Lowest order Hr^ (curl, f2)-conforming edge elements are widely used for the 
finite element Galerkin discretization of variational problems like (jl.ip . Then, for a 
solution u e {H^{n))^ with curlu G (iJ^(ri))'^ we can expect the optimal asymptotic 
convergence rate 

(1-3) \\^-^h\\„(^curin)<CN;;'^\ 
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jy(curl, il) : Sobolev space of square integrable vector fields on C M^^ with 
square integrable curl 

i?ri3 (curl, i7) : vector fields in i?(curl, £7) with vanishing tangential compo- 
nents on r^) C dil 

A4, T : tetrahedral finite element meshes, may contain hanging nodes 

A/'(7W) : set of vertices (nodes) of a mesh M 

S{M) : set of edges of a mesh M 

Pk, Pm '■ shape regularity measures 

h : - local meshwidth function for a finite element mesh 

- (as subscript) tag for finite element functions 
U(A^) : lowest order edge element space on A4 

hE ■ nodal basis function of U(A^) associated with edge E 

V{M) : space of continuous piecewise linear functions on M 

V2{M.) : quadratic Lagrangian finite element space on A4 

V2{M) : quadratic surplus space, see (|2.19p 

bp : nodal basis function of V{M) ("tent function") associated with 

vertex p 

25x(A^) : set of nodal basis functions for finite element space X on mesh 
M 

TLfi ■ nodal edge interpolation operator onto XJ{M), see ()2.7|) 

Th ■ vertex based piecewise linar interpolation onto V{M) 

Pp : space of 3-variate polynomials of total degree < p 

U(A^), V{A4): finite element spaces oblivious of zero boundary conditions 

^ : nesting of finite element meshes 

i{K) : level of element K in hierarchy of refined meshes 

u!i : refinement zone, see (|4.ip 

S; : refinement strip, see (I5.57P 

58y, 53u : sets of basis functions supported inside refinement zones, see 
611) 

Qh : quasi-interpolation operator, Def. 15.11 

Table 0.1 
Important notation used in this paper 



on families of finite element meshes arising from global refinement. Here, u/j is the 
finite element solution, Nh the dimension of the finite element space, and C > 
does not depend on Nh- However, often u will fail to possess the required regularity 
due to singularities arising at edges/corners of dfl and material interfaces [22,23]. 
Fortunately, it seems to be possible to retain l|1.3p by the use of adaptive local mesh 
refinement based on a posteriori error estimates, see [10, 55] for theory in iJ^-setting, 
[7,17] for numerical evidence in the case of edge element discretization, and [8,34,52] 
for related theoretical investigations. 

We also need ways to compute the asymptotically optimal finite element solution 
with optimal computational effort, that is, with a number of operations proportional 
to Nh- This can only be achieved by means of iterative solvers, whose convergence 
remains fast regardless of the depth of refinement. Multigrid methods are the most 
prominent class of iterative solvers that achieve this goal. By now, geometric multigrid 
methods for discrete if (curl, r2)-elliptic variational problems like (|l.ip have become 
well established [20,29,54,58]. Their asymptotic theory on sequencies of regularly 
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refined meshes has also matured [2,25,29,31,51]. It confirms asymptotic optimality: 
the speed of convergence is uniformly fast regardless of the number of refinement 
levels involved. In addition, the costs of one step of the iteration scale linearly with 
the number of unknowns. 

Yet, the latter property is lost when the standard multigrid correction scheme is 
applied to meshes generated by pronounced local refinement. Optimal computational 
costs can only be maintained, if one adopts the local multigrid policy, which was 
pioneered by A. Brandt et al. in [5], see also [41]. Crudely speaking, its gist is to 
confine relaxations to "new" degrees of freedom located in zones where refinement 
has changed the mesh. Thus an exponential increase of computational costs with 
the number of refinement levels can be avoided: the total costs of a V-cycle remain 
proportional to the number of unknowns. An algorithm blending the local multigrid 
idea with the geometric multigrid correction scheme of [29] is described in [54]. On 
the other hand, a proof of uniform asymptotic convergence has remained elusive so 
far. It is the objective of this paper to provide it, see Theorem 14.21 

We recall the key insight that (jl.ip is one member of a family of variational 
problems. Its kin is obtained by replacing curl with grad or div, respectively. All 
these differential operators turn out to be incarnations of the fundamental exterior 
derivative of differential geometry, cf. [29, Sect. 2]. They are closely connected in the 
deRham complex [3] and, thus, it is hardly surprising that results about the related 
i7p^ (rj)-elliptic variational problem, which seeks u € i/p^ (f?) such that 

(1.4) (gradw,gradw)^2(o) + (w,w)L2(f^) = (/,'y)i2(o) e i/^t^ (^2) , 

prove instrumental in the multigrid analysis for discretized versions of p.ip . Here 
Hp^{n) is the subspace of H^{il) whose functions have vanishing traces on F^. 

Thus, when tackling p.ip . we take the cue from the local multigrid theory for 
(|1.4p discretized by means of linear continuous finite elements. This theory has been 
developed in various settings, cf. [5,11,14,15,62]. In [1] local refinement with hanging 
nodes is treated. Recently, H. Wu and Z. Chen [60] proved the uniform convergence of 
local multigrid V-cycles on adaptively refined meshes in two dimensions. Their mesh 
refinements are controlled by a posteriori error estimators and carried out according 
to the "newest vertex bisection" strategy introduced, independently, in [6,40]. 

As in the case of global multigrid, the essential new aspect of local multigrid 
theory for (|l.ip compared to (|1.4p is the need to deal with the kernel of the curl- 
operator, cf. [29, Sect. 3]. In this context, the availability of discrete scalar potential 
representations for irrotational edge element vector fields is pivotal. Therefore, we 
devote the entire Sect. [2] to the discussion of edge elements and their relationship 
with conventional Lagrangian finite elements. Meshes with hanging nodes will receive 
particular attention. Next, in Sect. [3] we present details about local mesh refinement, 
because some parts of the proofs rest on the subtleties of how elements are split. The 
following Sect, m introduces the local multigrid method from the abstract perspective 
of successive subspace correction. 

The proof of imiform convergence fTheorem l4.2p is tackled in Sects. [5] and [51 which 
form the core of the article. In particular, the investigation of the stability of the local 
multilevel splitting requires several steps, the first of which addresses the issue for the 
bilinear form from (II. 4p and linear finite elements. These results are already available 
in the literature, but are re-derived to make the presentation self-contained. This 
also applies to the continuous and discrete Helmholtz-type decompositions covered in 
Sect. 15.31 Many developments are rather technical and to aid the reader important 
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notations are listed in Table 10.11 Eventually, in Sect. [71 we report two numerical 
experiments to show the competitive performance of the local multigrid method and 
the relevance of the convergence theory. 

Remark 1.1. In this article we forgo generality and do not discuss the more 
general bi-linear form 

(1.5) a(u,v) (acurlu,curlv)^2(o) + (/3u, v)^2(j^) , Vu, v e iJr^, (curl, fi) , 

with uniformly positive coefficient functions a,f3£ L°°(f2). We do this partly for the 
sake of lucidity and partly, because the current theory cannot provide estimates that 
are robust with respect to large variations of a and (3, cf. [33]. We refer to [63] for 
further information and references. 

2. Finite element spaces. Whenever we refer to a finite element mesh in this 
article, we have in mind a tetrahedral triangulation of O, see [19, Ch. 3]. In certain 
settings, it may feature hanging nodes, that is, the face of one tetrahedron can coincide 
with the union of faces of other tetrahedra. Further, the mesh is supposed to resolve 
the Dirichlet boundary in the sense that Fd is the union of faces of tetrahedra. The 
symbol A4 with optional subscripts is reserved for finite element meshes and the sets 
of their elements alike. 

We write h e L°°{n) for the piecewise constant function, which assumes value 
hj( := diam(if) in each element K E A4. The ratio of diam(i^) to the radius of the 
largest ball contained in K is called the shape regularity measure pk [19, Ch. 3, §3.1]. 
The shape regularity measure pM of A4 is the maximum of all pK, K G M. 

2.1. Conforming meshes. Provisionally, we consider only finite element 
meshes Ai that are conforming, that is, each face of a tetrahedron is either con- 
tained in Oft or a face of another tetrahedron, see [19, Ch. 2, § 2.2]. In particular, 
this rules out hanging nodes. Following [12,43], we introduce the space of lowest order 
i?rj3 (curl, r2)-conforming edge finite elements, also known as Whitney- 1-forms [59], 

(2.1) U(7W) := {v,, e i?rB(curl,r2) : y K e M : 3 a,b e : 

(2.2) -Vhix) = a + bx X, x e K} . 

For a detailed derivation and description please consult [30, Sect. 3] or the monographs 
[13,42]. Notice that curlU(A^) is a space of piecewise constant vector fields. We also 
remark that appropriate global degrees of freedom (d.o.f.) for \J{A4) are given by 

(2.3) ^ f , EeSiM), 
{ ^ JE^h-ds ' 

where S{M) is the set of edges of M not contained in F^. We write *Bu(A^) for the 
nodal basis of U(A^) dual to the global d.o.f. (12. 3|) . Basis functions are associated 
with active edges. Hence, we can write ^u{M) — {hE}EeSiM)- The support of the 
basis function is the union of tetrahedra sharing the edge E. We recall the simple 
formula for local shape functions 

(2.4) hE\K = K grad - grad A,; E = [a^, aj] C K 



for any tetrahedron K E Ai with vertices a^, i = 1, 2, 3, 4, and associated barycentric 
coordinate functions A^. 
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The edge element space 'U{A4) with basis ^u{Ai) is perfectly suited for the finite 
element Galerkin discretization of (|l.ip . The discrete problem based on U(7W) reads: 
seek u;i G \J{A4) such that 

(2.5) (curluft,curlv,0^2(o) + i^h,Vh)L2(n) = (f'V/0L2(j2) Vv,i G V{M) . 

The properties of \J{M) will be key to constructing and analyzing the local multigrid 
method for the large sparse linear system of equations resulting from (|2.5p . Next, we 
collect important facts. 

The basis ^u{M) enjoys uniform L^-stability, meaning the existence of a con- 
stanll3 C = C{pm) > such that for aU v,i J2 c^e'^e e U(X), G M, 

Ee£{M) 

(2.6) C-iK||^.(,,)< Yl ^l\\^E\\l2^n)<Chh\\l2^n) ■ 

Ee£{M) 

The global d.o.f. induce a nodal edge interpolation operator 
r dom(n,0 C i?ro(curl,17) ^ IJ{M) 

(2-7) V ^ E (/^v.d.-').b^,. 

Obviously, Tlh provides a local projection, but it turns out to be unbounded even 
on {H^{fl))^. Only for vector fields with discrete rotation the following interpolation 
error estimate is available, see [30, Lemma 4.6]: 

Lemma 2.1. The interpolation operator Tlh is bounded on G 
(iJ^(r2))'^, curl^f G curlU(7W)}c {H^{il))^, and for any conforming mesh there 
is C = C{pm) > such that 

||/t-i(/rf-nh.)*||^ <C|*|^i(f^) V* G (i^l(^^))^ curl* g curlU(7W) . 



If n is homeomorphic to a ball, then grad 71^(51) — i?(curlO, £7) := {v G 
/f(curl, fi), curlv — 0}, that is, H^{n) provides scalar potentials for J?(curl, 17). 

To state a discrete analogue of this relationship we need the Lagrangian finite ele- 
ment space of piecewise linear continuous functions on A4 

(2.8) ViM) := {uh G H^^in) : Uh\K e Fi{K) e M} , 

where Pp{K) is the space of 3-variate polynomials of degree < p on K. The global 
degrees of freedom for V{A4) boil down to point evaluations at the vertices of Ai 
away from Td (set Af{Ai)). The dual basis of "tent functions" will be denoted by 
^v{M) = {bp} 

peUiM)- unconditional L^-stability is well known: with a universal 

constant C > we have for all ut — ctp^p G y{M)^ ctp G K, 

peAf{M) 

(2.9) < Yl 4\\bp\\Un)<C\\uH\\l^n) ■ 

peM(M) 



^The symbol C will stand for generic positive constants throughout this article. Its value may 
vary between different occurrences. We will always specify on which quantities these constants may 
depend. 
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For the nodal interpolation operator related to *By we write 2/,, : dom(2/i) C 
_ffp^ (il) 1-^ V{Ai). Recall the standard estimate for linear interpolation on conforming 
meshes (i.e., no hanging nodes allowed), [19, Thm. 3.2.1], that asserts the existence 
of C = C{k,pM) > such that 

(2.10) \\h''-^ild -Ih)u\\j^,^^^^ < C\u\H2^n) yueH\n)nH^jn), fee {0,1}. 

Obviously, gradP^(AI) C U(A^), and immediate from Stokes theorem is the crucial 
commuting diagram property 

(2.11) Till o grad = gradol/j on dom(X/i) . 

This enables us to give an elementary proof of Lemma 1 2. II 

Proof, [of Lemma 12. 1| Pick one K Ai and, without loss of generality, assume 
^ K. Then define the lifting operator, cf. the "Koszul lifting" [3, Sect. 3.2], 

(2.12) w ^ Cw , C-w{x) |w(a;) x x , x £ K . 
Elementary calculations reveal that for any constant vectorfield w e {Vo{K))^ 

(2.13) curl£w=:w, 

(2.14) ll^w||^,(^) </ik||w||^,(^) , 

(2.15) CyveV{K). 

The continuity (|2.14p permits us to extend C to {L'^{K))^. 

Given * e {H^{K))^ with curl* = const^, by (|2.15p we know £curl* G 
{Vi{K))'^. Thus, an inverse inequality leads to 

(2.16) [/: curl <C/i^i II/: curl *||^.(^) < C ||curl *||^,(^) , 
with C = C{pk) > 0. Next, ((2TT3)) imphes 

(2.17) curl(* - Z:curl*) = 3peH^{K): * - £ curl * = gradp . 

From (|2.16p we conclude that p G H^{K) and |p|^2(^) < C|*|^i^^-). Moreover, 
thanks to the commuting diagram property we have 

(2.18) ^ -Uh^ curl ^ ~UhC curl * + grad(p - Ihp) , 

=0 by 1 I2T5I 

which means, by the standard estimate (|2.10|) for linear interpolation on K, 

II* -n?,*||i2(^) = \p~Ihp\Hi(K) < ChK\p\H2(K) < '^^K\^\m(K) ■ 

Summation over all elements finishes the proof. □ 

As theoretical tools we need "higher order" counterparts of the above finite ele- 
ment spaces. We recall the quadratic Lagrangian finite element space 

(2.19) V2{M) := {uh e H^^iQ) : Uh\K £ P2(i^) WK e M} , 
and its subspace of quadratic surpluses 

(2.20) V2{M) := {uh e V2{M) : IhUh = 0} . 
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This implies a direct splitting 

(2.21) V2{M) =V{M)®V2{M) , 

which is unconditionally _ff "'^-stable: there is a C = C{pm) > such that 

(2.22) C'^V/i|?/i(o) ^ - Ih)uh\%i^^-^ + llhUhln^^^ <C\uh\%i^n-^ , 
for all u/i e V2{M). 

Next, we examine the space {V{A4))^ of continuous piecewise linear vector fields 
that vanish on r^j. Standard afhnc equivalence techniques for edge elements, see [30, 
Sect. 3.6], confirm 

(2.23) 3C^C{pm)>0: ||n,*,||^.(,,) < C ||vI/,|1^,(,,) y^,,e{V{M) f. 

Lemma 2.2. For all yf^ e {V{M))^ we can find Vh G V2{M) such that 

"^h = Uh'^h +gra.dvh , 

and, with C = C{pm) > 0, 

C'' Il*/^lli2(n) < l|n,.*ft||i2(o) + ||gradiJft.||^.(f,) < • 

For the proof we rely on a very useful insight, which relieves us from all worries 
concerning the topology of Q: 

Lemma 2.3. //v e ifpi, (curl 0,17) and Il/iV = 0, then v e gradi/p^ (17). 

Proof. Since the mesh covers fl, the relative homology group iJi(i7;r£)) is gen- 
erated by a set of edge paths. By definition (|2.3p of the d.o.f. of U(A^), the path 
integrals of v along all these paths vanish. As an irrotational vector field with van- 
ishing circulation along a complete set of Fj^-relative fundamental cycles, v must be 
a gradient. □ 

Proof, [of Lemma 12. 2| Given Vf/j G {V{Ai))^, we decompose it according to 

(2.24) = Uh^h + {Id - n,0*ft . 

— :grad v/^ 

Note that curl{Id — H/J'S'/i is piecewise constant with vanishing flux through all 
triangular faces of Ai. Then Stokes' theorem teaches that curl{Id — Ilh)^h = 0. 

By the projector property of Tlh, {Id — n/,,)^/i satisfies the assumptions of 
Lemma 12.31 Taking into account that, moreover, the field is piecewise linear, it is 
clear that {Id — Y\.h)^h = grad?/; with 'tp G V2{M.). Along an arbitrary edge path 7 
in M. we have J^{Id — II/J'S'/i • ds = so that ip attains the same value (w.l.o.g. = 0) 
on all vertices of A4. The stability of the splitting is a consequence of (|2.23p . □ 

By definition, the spaces U(A^) and V{Ai) accommodate the homogeneous 
boundary conditions on F^). Later, we will also need finite element spaces oblivi- 
ous of boundary conditions, that is, for the case Fu = 0. These will be tagged by a 
bar on top, e.g., U(A^), V{M), etc. The same convention will be employed for notions 
and operators associated with finite element spaces: if they refer to the particular case 
Td = 0, they will be endowed with an overbar, e.g. Ilh, 1h, Su(A^), A/'(7W), etc. 
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2.2. Meshes with hanging nodes. Now, general tetrahedral meshes with 
hanging nodes are admitted. We simply retain the definitions (|2.8p and (|2.19p of the 
spaces V{A4) and V2{Ai) of continuous finite element functions. Degrees of freedom 
for V{A4) are point evaluations at active vertices of A^. A vertex is called active, if 
it is not located in the interior of an edge/face of or on F/j. A 2E0 illustration is 
given in Fig. 12.11 

I I 



A 



Mo Ml M2 Ms 

Fig. 2.1. Active vertices (*) of 2D triangular meshes with hanging nodes, Q =]0, F^i = dQ. 
In A^i,jVf2, AI3 active edges are marked with green arrows. 




The values of a finite element function at the remaining ("slave") vertices are 
determined by recursive afhne interpolation. A dual nodal basis *By(A^) and corre- 
sponding interpolation operator can be defined as above. 

In principle, the definition (|2.ip of the edge element space could be retained on 
non-conforming meshes, as well. Yet, for this choice an edge interpolation operator 
Tlh that satisfies the commuting diagram property (|2.1ip is not available. Thus, we 
construct basis functions directly and rely on the notion of active edges, see Fig. 12.11 

Definition 2.4. An edge of M is active, if it is an edge of some K G M, not 
contained inTu, and connects two vertices that are either active or located onTo- 

We keep the symbol S{M) to designate the set of active edges of A^. To each 
E G £{M) we associate a basis function b^, which, locally on the tetrahedra of M, is 
a polynomial of the form (|2.2p . In order to fix this basis function completely, it suffices 
to speficify its path integrals (|2.3p along all edges of M. In the spirit of duality, we 
demand 



(2.25) / hs-ds^ 



' 1 ,ifF^E, 



F ' \0 ,ifFe£{M)\{E} 



For the non-active ( "slave" ) edges of M the path integrals of b^; (subsequently called 
"weights") are chosen to fit ()2.1ip . keeping in mind that ^\j{M) := {i>E}Ee£(M)^ 
and that the d.o.f. and Ilh are still defined according to (|2.3p and l|2.7p . respectively. 
Uhimately, we set U(X) Span {Q3u(A^)}- 

Let us explain the policy for setting the weights in the case of the subdivided 
tetrahedron of Fig. 12. 21 with hanging nodes at the midpoints of edges, which will turn 
out to be the only relevant situation, cf. Sect. 15.31 Weights have to be assigned to the 
"small edges" of the refined tetrahedron, some of which will be active, and some of 
which will have "slave" status, see the caption of Fig. 12.21 



^For ease of visualization, we will often elucidate geometric concepts in two-dimensional settings. 
Their underlying ideas are the same in 2D and 3D. 
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We write the direction vectors of slave edges as linear combinations of active 
edges, for instance, 

9i - 92 = |(P4 +P3) - 5(^2+^3) = 5(^4-^2) , 

P5 - 94 = P5 - liPl + P3) ^P^-Pi + \{Pi- P3) ' 

94 - 93 = Upi+Ps) - 5(P4 +P2) ^ -P2) + 5(^3-^4) ■ 

In a sence, we express slave edges as "linear combinations" of active edges. In a 
different context, this policy is explained in more detail in [26]. 




Fig. 2.2. Subdivided tetrahedron, active vertices (•) P]^,...,pg, slave vertices (o) qi,...,q^, 
active edges [pi.Pg], [Pi.Pe], [Pi^Ps]' fe.Pel' fe-Pal' [P2.P4]. [P3'P4]. [PS'Pe]- [P1.P3]. '^^"■'"^ 
edges [pi, 94], [q4,,Ps], [P2,q2]. [92>P3]. [P3'9i], [9i,P4]. [Pe-^J. ['J2:q4]> [P6:92]. [9i.92]. [92.93]. 
[91.93]' [91.P5]. [93. P5]. [94. P5] 

The coefficients in the combinations tell us the weights. For example, for the active 
edge E = [^3,^4] in Fig. 12.21 thev are given in Table [27T] Using these weights and the 
formula (|2.4p . b^; can be assembled on the tetrahedron by imposing (see Table \TJ\ for 
notations) 



be; • ds 



for any contributing slave edge S , 
for all other (slave) edges. 



S £ { "small edges" } 



Firstly, the procedure for the selection of weight guarantees that gradF(7W) C 
U(A^). For illustration, we single out the gradient of the nodal basis func- 
tion belonging to vertex pg in Fig. 12.21 Its path integral equals 1 along the (ori- 
ented) edges [pi,P5], [P3,P5], [pg,P5], [94,^5], [93,^5], [^1,^5], and vanishes on all 



10 



R. Hiptmair and W.-Y. Zheng 



Slave edge S 


[qi,P4] 




[92,93] 


[94:^5] 


[93>94] 


weight ws 


1 

2 


1 

2 


1 
2 


i 
2 


1 

2 



Table 2.1 

VKeij/iis /or slave edges in Fig. |g. 2l relative to active edge E = [p3,p4]. Only slave edges with 
non-zero weights are listed. 



other edges. Hence we expect 

(2.26) w,, = b[p^^p^] + b[p^^^p^] + b[p^,p^] . 

This can be verified through showing equahty of path integrals along slave edges. 
We take a close look at the slave edge [94,^5]. By construction the basis functions 
belonging to active edges satisfy 

J b[p^^p^] •ds = 1 , J b[p^^pj •ds = 0, J b[p^,pj -ds = , 

[<li,P5] lli^Pr,] lli^Ps] 

J b[p,,p3] -ds^ -i , J b[p^^p^] -ds^O J b[p^^p^] •ds = 0, 



lli.Ps] lli^PrJ lli-^P^] 



194. P5I 

Then, evidently, 
1 = y w/i • ds 



'^[Pi.Ps] ■ d'?+ J ''[Ps.Psl • ds + J b[pg^p^] • ds = 1 + + . 

['J4.Pb] [q4.P5] ['J4'P5] 

The same considerations apply to all other slave edges and (|2.26p is established. 
Secondly, the construction ensures the commuting diagram property (|2.1ip : again 
appealing to Fig. 12.21 we find, for example. 



J gradI,iU • ds^Ihuiq^) - Ihu{q^) = 
1Wp,) + „(P.,)-1Hp,) + „(p,)) = 1 / grad..d,r+l / .™d„.d,r. 

[P3.P4I [Pl^P2] 

i j grad u • ds + i j grad u ■ (is + j j grad u • ds . 

[Ps^Pd] [Pl.Ps] [P6.P2] 

In words, combining the path integrals of grad u along active edges with the relative 
weights of the slave edge [^3, yields the same result as evaluating the path integral 
of the gradient of the intcrpolant Ihu along [qig, q^. 

The definitions (j2.19|) and ()2.20|) also carry over to meshes with hanging nodes. 
This remains true for the splitting asserted in Lemma 2.2. However, though the al- 
gebraic relationships like (j2.1ip remain valid, the estimates and norm equivalences of 
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the previous section do not hold for general families of meshes with hanging nodes. 
This entails restrictions on the location of hanging nodes, whose discussion will be 
postponed until Sect. 13.21 cf. Assumption 13 . . ll 

Remark 2.5. Our presentation is confined to tetrahedral meshes and lowest order 
edge elements for the sake of simplicity. Extension of all results to hexahedral meshes 
and higher order edge elements is possible, but will be technical and tedious. 

3. Local mesh refinement. We study the case where the actual finite element 
mesh of 17 has been created by successive local refinement of a relatively uniform 
initial mesh A4q. Concerning Ad^ and A4q the following asumptions will be made: 

1. Given A4q and Mh we can construct a virtual refinement hierarchy of i + 1 
nestecjl tetrahedral meshes, L e N: 

(3.1) Mo Ml < M2 -< ■■■ -(ML=Mh. 

Please note that the virtual refinement hierarchy may be different from the 
actual sequence of meshes spawned during adaptive refinemen10. 

2. Inductively, we assign to each tetrahedron K e Mi a level £{K) e No by 
counting the number of subdivisions it took to generate it from an element 
oi Mq. 

3. For all < / < L the mesh Mi+i is created by subdividing some or all of the 
tetrahedra in {K e Mi : i{K) = I}. 

4. The shape regularity measures of the meshes Mi are uniformly bounded 
independently of L. 

Refinement may be local, but it must be regular in the following sense, cf. [46, 
Sect. 4.2.2] and [60]: we can find a second sequence of nested tetrahedral meshes of 

(3.2) Mo^Mq -< Ml -< M2 -< ■■■ <Ml ■ 

that satisfies 

1. Ml -< Ml and {K e Mi : £{K) = 1} C Mi, I = 0, . . . , L, 

2. that the shape regularity measure pj^^ is bounded independently of I, 

3. and that there exist two constants C > and < 9 < 1 independent of / and 
L such that 

(3.3) C-^e*' <hK <C9^ VK eMi, 0<1<L. 

This means that the family {Mi}i is quasi-uniform. Hence, it makes sense to 
refer to a mesh width hi := maxj/i/j-. A' G A^;} of A^;. It decreases geomet- 
rically for growing /. 
Our analysis targets two popular tetrahedral refinement schemes that generate 
sequences of meshes that meet the above requirements. 

3.1. Local regular refinement. This scheme produces Mi+i by splitting some 
of the tetrahedra of the current mesh Mi into eight smaller ones, possibly creating 
hanging nodes in the process [1]. An illustrative 2D example with hanging nodes is 
depicted in Figure [57T] The accompanying sequence {Mi}a<i<L is produced by global 
regular refinement, which implies (j3.3p with 9 = ^. Uniform shape- regularity can also 
be guaranteed for repeated regular refinement of tetrahedra, see [9]. 



^ two finite element meshes and T arc nested, M -< T, if every element of M is the union of 
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Mq Ml M2 M3 = Mh 

Fig. 3.1. Virtual refinement hierarchy for 2D triangular meshes. The quasi-uniform sequence 
{•^!}o<i<L sketched in blue. Elements of Mi eligible for further subdivision are marked yellow. 



The meshes occurring in the virtual refinement hierarchy need not agree with the 
meshes that arise during adaptive refinement in an actual computation. Yet, given 
Mh, the virtual refinement hierarchy can always be found a posteriori. Write Mhicr 
for the union of all tetrahedra ever created during the refinement process. Then, for 
<l < L, define 

,^ , . ,, r„ ,, £(K) < I and if does not contain a 1 

(3-4) ■= ^ -^^^ ■ K' e X^M.. \ {K} with i{K') <l ]■ 

Using the construction of finite element spaces detailed in Sect. I2.2[ the local 
multigrid algorithms can handle any kind of local regular refinement. Yet, convergence 
may degrade unless we curb extreme jumps of local meshwidth. Thus, we assume the 
following throughout the remainder of this paper. 

Assumption 3.0.1. Any edge oj Mh may contain at most one hanging node. 

This will automatically be satisfied for all meshes Mi of the virtual refinement 
hierarchy. Consequently, hanging nodes can occur only in a few geometric configu- 
rations, one of which is depicted in Fig. 12.21 This paves the way for using mapping 
techniques and scaling arguments, see [30, Sect. 3.6], which confirm the following gen- 
eralization of results of Sect. 12.11 Of course, we rely on the constructions of finite 
element spaces and interpolation operators described in Sect. 12.21 

Proposition 3.1. Under Assumvtion \3.0.1\ the L^-stahility of bases, see (|2.6p . 
(|2.9[) . carries over uniformly to meshes created by local regular refinement. So do 
Lemmas\2l^\2M and Estimates ^(TF^ . jMS]). 

Summing up. Assumption 13.0.11 makes it possible to use the results obtained in 
Sect. 12.11 in the case of local regular refinement as well. To avoid a proliferation of 
labels, we arc going to quote the statements from Sect. I2.1l even when we mean their 
generalization to meshes with hanging nodes. 

3.2. Recursive bisection refinement. This procedure involves splitting a 
tetrahedron into two by promoting the midpoint of the so-called refinement edge 
to a new vertex. Variants of bisection differ by the selection of refinement edges: The 
iterative bisection strategy by Bansch [4, 6] needs the intermediate handling of hang- 
ing nodes. The recursive bisection strategies of [36,38,57] do not create such hanging 



elements of T. 

*For the local multigrid algorithm examined in this article the implementation must provide 
access to the virtual refinement hierarchy. This entails suitable bookkeeping data structures, which 
are available in the ALBERTA package used for the numerical experiments in Sect. [7] 
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nodes and, therefore, are easier to implement. But for special Mq, the two recursive 
algorithms result in exactly the same tetrahedral meshes as the iterative algorithm. 
Since our implementation relies on the bisection algorithm of [36], we outline its bi- 
section policy in the following. For more information on bisection algorithms, we refer 
to [49,56]. 

For the recursive bisection algorithm of [36] , the bisections of tetrahedra are to- 
tally determined by the local vertex numbering of M.o, plus a prescribed type for 
every element in A^o- Each tetrahedron K is endowed with the local indices 0, 1, 2, 
and 3 for its vertices. The refinement edge of each element is always set to be the edge 
connecting vertex and vertex 1. After bisection of K, the "child tetrahedron" of K 
which contains vertex of if is denoted by Child [0] and the other one is denoted by 
Child[l]. The types of Child[0] and Child[l] are defined by 

type(Child[0]) = type(Child[l]) = (type(if) + 1) mod 3. 

The new vertex at the midpoint of the refinement edge of K is always numbered 
by 3 in Child[0] and Child[l]. The four vertices of K are numbered in Child[0] and 
Child[l] as follows (see Fig.^: 

InChild[0]: (0,2,3)^(0,1,2), 

In Child[l] : (0, 2, 3) ^ (0, 2, 1) , if type(iC) = 0, 

In Child[l] : (0, 2, 3) ^ (0, 1, 2) , if type(X) > 0. 

This recursive bisection creates only a small number of similarity classes of tetrahedra, 
see [36,49,57]. 

Fig 13.31 shows a 2D example of the recursive bisection refinement (the algorithm 
for 2D case is called "the newest vertex bisection" in [41]). Similar to the 3D algorithm, 
for any element its three vertices are locally numbered by 0, 1, and 2, its refinement 
edge is the edge between vertex and 1. The newly created vertex in the two children 
of K are numbered by 2. In the child element containing vertex of K, vertex and 
2 oi K are renumbered by 1 and respectively. In the other child element, vertex 1 
and 2 oi K are renumbered by and 1 respectively. 

In order to keep the mesh conforming during refinements, the bisection of an edge 
is only allowed when such an edge is the refinement edge for all elements which share 
this edge. If a tetrahedron has to be refined, we have to loop around its refinement 
edge and collect all elements at this edge to create an refinement patch. Then this 
patch is refined by bisecting the common refinement edge. A more detailed discussion 
can be found in [36] . 

For any mesh A4i an associated "quasi- uniform" mesh A4i according to (|3.2p . 
M.I -< A4i, is obtained as follows: the elements in {K G A4i : £{K) < 1} undergo 
bisection until i{K) = I for any K £ 

We still have to make sure that the recursive bisection allows the definition of 
a virtual refinement hierarchy. Thus, let Aih — Ml generated from the initial 
mesh A4q by the bisection algorithm in [36]. Denote by A^hicr the set of all tetrahedra 
created during the bisection process, i.e., for any K e A^hicr, there is a if' G Mh 
such that either K' — K or K' is created by refining K. Then, the virtual meshes 
Ml, < I < L can again be defined according to ()3.4|) . 

In the following, we are going to prove that each A4i is a conforming mesh, that 
is, no hanging nodes occur in A4i, < I < L. The proof depends on some mild 
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Fig. 3.2. Bisection of tetrahedra in the course of recursive bisection. Assignment of types to 
children 



assumptions on Aio (see assumptions (Al) and (A2) in [36]), which wiU be taken for 
granted. 

Lemma 3.2. [36, Lemmas 2,3] Let T,T' G Ai^ be a pair of tetrahedra sharing a 
face F — K n K' . It holds true that 

1. if T contains the refinement edge of T' and vice versa, then they have the 
same refinement edge, 

2. if F contains the refinement edges of both K and K' , then i{K) — i{K'), 

3. if F contains the refinement edge of K , but does not contain the refinement 
edge of K' , then i{K) = i{K') + 1, 

4- if F does not contain the refinement edges of K and K' , then £{K) = £{K')- 
Lemma 3.3. The meshes Mi, < I < L, according to p.4p are conforming 
meshes. 

Proof. We are going to prove the lemma by backward induction starting from 
I = L. Since Ml = Mh is conforming, for any K e Ml satisfying t{K) = L, there 
exists a brother of K, denoted by K' G Ml, such that i{K') = L and Kp := KUK' G 
Ml-1- Here Kp is called the parent of K and K' with £{Kp) = i - 1 (see Fig. 13. 4p . 

Let E be the refinement edge of Kp. By the recursive bisection algorithm, E must 
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M4 




6 



Fig. 3.3. Virtual refinement hierarchy for 2D triangular meshes emerging in the course of 
successive local newest vertex bisection refinement of Mo from Fia. \3. ll Accompanying quasi-uniform 
meshes outlined in blue, maximally refined triangles marked yellow. 



be the common refinement edge of all tetrahedra in the refinement patch: 

By Lemma [3. 2 [ i{Kp) — L — 1 ioi any C Pe and the midpoint of E, denoted by 
^now, is the unique new vertex of A^l in Pe- We conclude that 

Pe = [J{^ : K E Ml, i{K) = L, and ^ncw is a vertex of K }. 

Coarsen the sub-mesh A^£|p^ by removing the vertex A^cw and all edges related to it 
and adding E to this patch. Thus a conforming sub-mesh A^l-1|p^ is obtained. Do 
the above coarsening process for every element K e Ml with i{K) = L. This proves 
that Ml-1 is conforming. 

Finally, an induction argument confirms that Mi is conforming, I = L — 2, ■ ■ ■ ,1. 

□ 

4. Local multigrid. To begin with, wc introduce nested refinement zones as 
open subsets of Q: 

(4.1) uji := inferiority {7? : K G Mh, 1{K) > O) C f7 , 

see Fig. 14.11 and Fig. 14.21 The notion of refinement zones allows a concise definition of 
the local multilevel decompositions of the finite element spaces V{Mh) and \J{Mh) 
that underly the local multigrid method. 



16 



R. Hiptmair and W.-Y. Zheng 




We introduce local multigrid from the perspective of multilevel successive sub- 
space correction (SSC) [61,62,64]. First, we give an abstract description for a linear 
variational problem 

(4.2) u e H : a(u, v) = f{v) Vw G , 

involving a positive definite bilinear form a on a Hilbert space H. The method is 
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Fig. 4.2. Refinement zones for the 2D refinement hierarchy of Figure [3. 3\ 



completely defined after we have provided a finite subspace decomposition 
,7 

(4.3) H ^"^Hj , Hj c H closed subspaces, j = 0, . . . , J, J e N . 

j=o 

Then the correction scheme implementation of one step of SSC acting on the iterate 

for m = 1, 2, • • • 

for j = 0,l,-- - 

Let Cj G Hj solve 

a (ej , ) = / {vj ) - a (u™"/ ,Vj) V i;j G iJ^ 

endfor 
endfor 

This amounts to a stationary linear iterative method with error propagation op- 
erator 

(4.4) E={I- Pj)il - Pj^i) •••(/- Po) , 

where Pj : H Hj stands for the Galerkin projection defined through 

(4.5) a{PjV,Vj) = a{v,Vj) yv^eHj. 

The convergence theory of SSC for an inner product a and induced energy norm 
ll'll^ rests on two assumptions. The first one concerns the stability of the space de- 
composition. We assume that there exists a constant Cstab independent of J such 
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that 



(4.6) inf{^||«,||^: ^z;, <Cstab|k||^ V« £ 

The second assumption is a strengthened Cauchy-Schwartz inequality, namely, there 
exist two constants < q < I and Corth independent of j and k such that 

(4.7) a{vj,vk) < Corth'?'''"^'' \\vj\\a IMa ^ Hj, Vk e Hk . 

The above inequality states a kind of quasi-orthogonality between the subspaces. 
From [61, Theorem 4.4] and [66, Theorem 5.1] we cite the following central convergence 
theorem; 

Theorem 4.1. Provided that (14. 6p and (j4.7p hold, the convergence rate of Algo- 
rithm SSC is bounded by 

(4.8) \\E\\1<1~- ^i— — wtth e = Corth^, 

'-'stab I -|- ^ ~ Q 

where the operator norm is defined by 

||£|U:^ .,u„ leil. 

veH,v^O \\V\\a 

The bottom line is that the subspace splitting (|4.3p already provides a full de- 
scription of the method. Showing that both constants Cstab from (j4.6p and Corth from 
(j4.7p can be chosen independently of the number L of refinement levels is the challenge 
in asymptotic multigrid analysis. 

In concrete terms, the role of the linear variational problem ()4.2p is played by 
(jl.ip considered on the edge element space U(7W/j), which replaces the Hilbert space 
H. To define the local multilevel decomposition of \J{Aih), we define "sets of new 
basis functions" on the various refinement levels 

^l:^^v{Mo), ^\,:={bhe'Bv{Mi): snppbhCcJi} , ^ ^ ^ ^ 

<8?j *Bu(A^o), ^[j -.^{hne^u (Ml) : supphhCuJi} , " 

A 2D drawing of the sets Q5y is given in Fig. 14.31 where = dil. Note that we 
also have to deal with V{A4h.), because, as suggested by the reasoning in [29], a local 
multilevel decomposition of XJ{A4h) has to incorporate an appropriate local multilevel 
decomposition of V{Aih)- 

Then, a possible local multigrid iteration for the linear system of equations arising 
from a finite element Galerkin discretization of a H^^ (r2)-elliptic variational problem 
boils down to a successive subspace correction method based on the local multilevel 
decomposition 

L 

(4.10) V{Mh) = V{Mo) + Y, Span {64 . 

Similarly, the local multilevel splitting of XJ{Mh) is based on the multilevel de- 
composition 

L L 

(4.11) V{Mh)=V{Mo)+Y^ Span{grad64 + ^ ^ Span{b4 . 
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These splittings induce SSC iterations that can be implemented as non-symmetric 
multigrid V-cycles with only one (hybrid) Gauss-Seidel post-smoothing step, see [29, 
Sect. 6]. Duplicating components of (|4.1ip results in more general multigrid cycles 
with various numbers of pre- and post-smoothing steps. 

The splitting (14. lip is motivated both by the design of multigrid methods for 
(|l.ip and U(A^) in the case of uniform refinement and local multigrid approaches 
to Hp^ (r2)-elliptic variational problems after discretization by means of linear finite 
elements [41,60]. The occurrence of gradients of "tent functions" bh in (|4.1ip is related 
to the hybrid local relaxation^ which is essential for the performance of multigrid 
in ff(curl, f2), see [29] for a rationale. A rigorous justification will emerge during 
the theoretical analysis in the following sections. It will establish the following main 
theorem. 

Theorem 4.2 (Asymptotic convergence of local multigrid for edge elements). 
Under the assumptions on the meshes made above and allowing at most one hanging 
node per edge, the decomposition ()4.1ip leads to an SSC iteration whose convergence 
rate is bounded away from 1 uniformly in the number L of refinement steps. 

5. Stability. First we tackle the stability estimate (|4.6p for the local multilevel 
decomposition ()4.10p . which is implicitly contained in ()4.1ip . 

5.1. Local quasi-interpolation onto V{A4). Quasi-interpolation operators 
are projectors onto finite element spaces that have been devised to accommodate two 
confiicting goals: locality and boundedness in weak norms [18,21,46,50,53]. As key 
tool they will be used in Sect. 15.21 and the proof of Lemma [5. 101 As in [46, Sect. 2.1.1], 
we resort to a construction employing local linear L^-dual basis functions. We follow 
the analysis of [50] that permits us to take into account Dirichlet boundary conditions. 

For a generic tetrahedron K define tpf ^ J = lj2,3,4, by L'^ {K)-duality to the 
barycentric coordinate functions A^, i = 1, 2, 3, 4, of K: 

(5.1) VfGFi(if): [ ^pf{x)X,{x)dx = d,j , z,je{l,...,4}. 

Jk 

Computing an explicit representation of the i/jj^ we find 

(5.2) C-'<\K\\\^f\\l^^^<C , C-i<|K||^,(^^<C, 

with an absolute constant C > 0. We can regard ip^^ as belonging to the j-th vertex 
of K. Thus, we will also write ipp , p G M{K), N{K) the set of vetices of K. 
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Let Ai be one of the tetrahedral meshes A4i ot Mi oi il. In order to introduce 
quasi-interpolation operators we take for granted some "node— ^ceh" -assignment, a 
mapping 77{M) ^ M, p ^ 7J{M) ^ Kp £ M. 

Definition 5.1. Writing {^p}pgjv-(A4) ^v{M), define the local quasi- 
interpolation operator 



(5.3) 



Qh 



L^{n) V{M) 



T.peM{M) Ik^ V-p " ix)u{x) dx ■ h. 



Analoguously, we introduce the local quasi-interpolation Q/j : L^(ri) i— > V{M). 

We point out that Qh respects u — Q onT because the sum does not cover basis 
functions attached to vertices on To- From (|5.ip it is also evident that both Qh and 
Qh are projections, for instance, 

(5.4) QhUh = Uh Wuh e V{M) . 

Moreover, they satisfy the foUowing strong continuity and approximation properties: 



Lemma 5.2. The quasi-interpolation operators from Def. I5.il allow the estimates 
(set To^iforQM) 



|Q/»m|Il2(o) < C'||w||i2(o) 



(5.5) 

3C = CipM) ■■ 
(5.6) 

(5.7) 

3C = C{pM,k): \\h-''{u-Qhu)\\^ <C\u\^,,^. V u £ h'' {n) n jn) 



and k — 1.2. 



Proof. [Part I] Continuity in L'^{n) is a simple consequence of the stability 
of the nodal bases Sy(A4) and of the Cauchy-Schwarz inequality: 



\Qhu 



iL^n)<C E \Qhuip)\' \\bp 
pei^{M) 



2 



ipp ^{x)u{x) dx 



peAfiM) 



with C = C{pm) > 0, because 



< C, too. □ 



The following estimate is instrumental in establishing continuity of Qh in (fl): 
Theorem 5.3 (Generahzed Hardy inquahty). 



3C ^Cin,TD)>0: J 



dist(a;, r^i) 
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Proof. By density it suffices to consider u e C°°{il), supp(w) n Td — 0- Using a 
partition of unity, we can confine the estimate to neighborhoods of To , in which 90 is 
the graph of a Lipschitz-continuous function. Thus, after bi-Lipschitz transformations, 
we need only investigate three canonical situations, see Fig. 15.11 

1. Td = {z = 0}, for which the ID Hardy inequality gives the estimate, see the 
proof of Thm. 1.4.4.4 in [27]. 

2. Td = {z = Q a X > 0}, which can be treated using polar coordinates in the 
{x, z)-plane and then integrating in y-direction: 



u{r, ip) 



dip r dr < 



TT du 

-jr[r, p) 

r oip 



d^.d.<.^/|grad,,.Pdxd. 



z>0 



3. To = {z = OAx>OAy> 0}, for which we obtain a similar estimate using 
spherical coordinates. 
This ends the proof. □ 

z z z 






V 




/ 














Vd = {z 


= A a; > 0} To 



D 



V 



dn\rD 



Fig. 5.1. Canonical situations to be examined in the proof of Thm. [5731 

Proof, [of Lemma [5.21 part II] In order to tackle the _ff ^(il)-continuity of Qh, we 
use that gradV^(A^) C U(A^) along with the stability estimate (|2.6p 



(5.8) ||gradQ,«||2,^^^ <C ^ {Q^,u{p) ~ Q,,u{q)f \\hE\\l.^^^ , 

E=[p,q]e£{M) 

with the notation {^E}Ee£(M) ■~ Su(-M)- 

(i) for the case E — [p,q] £ £{A4), p,q ^ Tjj, we adapt arguments from 
For any u e H^^{n), by (|5.ip . we have the identity 



\{Qhu){p) ~ {Qhu)iq)\ 



il}p''{x)i)q''{y){u{x) -u{y))dydx 



■4!p''{x)i'q''{y) I ^vsLdu{y + t{x - y)) ■ {x - y)dLT dydx 



Then split the innermost integral and transform 



f{v + t{x - y))dT ^ f{y + T{x-y))dT + f{x + T{y-x))dT 
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We infer 

\iQhu){p) - iQhu)iq)\ 



1 

< J J J \i'p''{x)\\iljq''{y)\\gradu{y + r{x-y))\\x-y\dydxdT 

The transformation formula for integrals reveals 

[ f{x + T{y-x))dy = T-^ [ f{z)dz, K' := x + t{K - x) . 

JK JK' 



Appealing to the bounds for ||V'j^||^2(;f): ll^i^lLi(i<-)' ^ ^ ^^m (|5.2p . the Cauchy- 
Schwarz inequality yields 



(5.9) \(QHu){p)~{Qhu){q)\\\hE\\L.^^^ 

I|b£;||i2(j^) 



< C 



imxi{\Kq\^,\Kp\^] 



T-3/2dr. 



\m{{nE)) ■ 



<C=C(pm) 



Here (51s) stands for the convex hull of all tetrahedra adjacent to the edge E. 
(ii) now consider E — [p, q\ G E{M), p eTd- Then, for any u e _ffp^ (il) 



\{Qnu){q) - {Qhu){p)\' = \{QHu){q)\' = 



tpq'' {x)u{x) dx 



dist{x,rD)i'q''{x)- "^^^ 



<C|p--q| 



<- 



c 



\P-Q\ 



u{x) 



dist(a;, To) 



dist(a;, Fd) 



da; 



■ da; 



2 ^ C\p-q\ 



Ka 



u[x) 



dist(a;, Fd) 



da; 



with (different) constants C — C{pm) > 0. 

Combining (jS.Sp . (|5.9p . using the finite overlap property of M in the form 

3C = C(pAi) : tt{£^ e S{M) : x e (fls)} <C \/x e Q , 

and appealing to Thm. [531 confirm |Qhu|//i(f;) < C'I'"Ihi(o)- Observe that the Hardy 
inequality makes the constant depend on f2 and F^i in addition. 

The quasi-interpolation error estimate ()5.7|) results from scaling arguments. Pick 
K E M, u e H^{il.) n H^^{Q), and write Iku G Fi{K) for the linear interpolant of 
u on K . Thanks to the projection property, we deduce as in Part I of the proof that, 
with C = C{pm)i 

||(/rf- Q/i)w||^2(^) = \\{Id-Qh){u-lKu)\\^2^K-^ < C\\u-lKu\\^2(^aK) 
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Here, we wrote fix '■— [J{K' : K' n K ^ 0}, and the final estimate can be shown 
by a simple scaling argument, cf. (|2.10p . Estimate (15. 7p for fc = 1 follows by scaling 
arguments and interpolation between the Sobolev spaces H^{VIk) and L^{Q.k)- □ 

5.2. Multilevel splitting of V{Ml)- In this section we first revisit the well- 
known [45,61,67] uniform stability of multilevel splittings of _ff ^(r2)-conforming La- 
grangian finite element functions in the case of mesh hierarchies generated by uniform, 
i.e. non-local, regular refinement. 

We take for granted a virtual refinement hierarchy (|3.ip of tetrahedral meshes as 
introduced in Sect. [3] and its accompanying quasi-uniform family of meshes (|3.2p . 

Owing to the inf in (|4.6p . it is enough to find a concrete family of admissible 
"candidate" decompositions that enjoys the desired L-uniform stability. We aim for 
candidates that fit the locally refined mesh hierarchy. 

The principal idea, borrowed from [46, Sect. 4.2.2], is to use a sequence 
of quasi-interpolation operators Q; : L^(ri) ^ V{Xii) based on a judiciously 
chosen node-^element-assignments. For M.i we introduce a "coarsest neighbor 
node—^element-assignment" : First, for any p G J\f{Mi), I = 1, . . . , L, we pick K G A4i 
such that 

e{K) = mm{e{K) : p eK, K E Mi} . 

Secondly, we select a "coarsest neighbor" Kp G Mi among those elements of Mi that 
are contained in K. This defines a mapping Af{Mi) i-^ Mi, p i— > Kp. We write Q; : 
L'^{Q) 1-^ V{Mi) for the induced quasi-interpolation operator according to Def. 15.11 
Next, we examine the candidate multilevel splitting 

L 

(5.10) Uh^QoUh + J2iQi-Qi-i )uh , Uh G V{Ml) ■ 

1=1 

Lemma 5.4. There holds, with a constant C > depending only on fl and the 
uniform bound for the shape regularity measures pj^^ , ^ I ^ L, 

L 

(5.11) |Qo"/.|Jfi(n)+IIV'||(Q;-Q(-iK||L2(o) <C\uh\liin) Vu,, G I^(A^l) ■ 

1=1 

Proof. We take the cue from the elegant approach of Bornemann and Yserentant 
in [11], who discovered how to bring techniques of real interpolation theory of Sobolev 
spaces [37], [39, Appendix B] to bear on (|5.10p . The main tools are the so-called K- 
functionals given by 

Kit,uf :=^^mf^^^|llM-u;11^2(fj)+i^lwl^2(o)} , 

KR3{t,uf := ^^inf^^^ - ^[[^^(Ra) + |u'Ih2(r3)} • 

The estimates (|5.6p and (I5.7P of Lemma 15.21 create a link between the terms in (|5.1ip 
and K{t,u): owing to (|5.5p and (|5.7p there holds for any u G L^{n) 

\\(^i~'^i-M\mn) ^ -Q'-i)("-^)llL^(fi) + ll(Q'-Q'-i)^IL2(n) 

< C{\\u - -I- hf\w\„.^,,^) yw G H^n) . 
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Here and below the generic constants C may depend on shape regularity max pj^^ 
and the (quasi-uniformity) constants in (13. 3p . We conclude 

(5.12) \\(^i-Qi_^)u\\^^,^^^^<CK{hluf yuGL^m, 
which implies 

L L 

(5.13) |Qo?/|^,(,,)+5]v'l|(Q/-Qi-i)"llL^(ii) <c^{l^l?/ia^)+EV^('^''")} • 

1=1 1=1 

Let u G H^{M.'^) be the Sobolev extension of u such that, with C = C{il) > 0, 

2|n = u and |u|^i(r3) < C'|w|jji(j^). 
Define the Fourier Transform of u by 

By the equivalent definition of Sobolev-norms on 

|S|^,(R3)« / ||Pl52(OI'd| , z = 0,l, 

we have 

because the infimum is attained for [39, Thni. B7] w{^) — u{^)/{l + Since 



inf 

uG/i"2(R3) 



inf {||M-u'lli2(s:2) < KT^^^it^uf 



we deduce that 



(5.14) 5:V'i^(/^^«)2<C^ / «^|^2(0pd^ 

<csup{x: ,f''g|' } / i^n^i(ordi 

< C |u|//i(R3) < C'l'«l/i-i(r2), 

where we have used assumption p.3p . The proof is finished by combining (|5.13p and 
(I5T41) . □ 

Now we restrict ourselves to Uh G V^(A^ft,). Then, thanks to the particular design of 
the node^element-assignment underlying Q/, the terms in the decomposition (|5.10p 
turn out to be localized. 
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Lemma 5.5. For all Uh G V{Mh) and < I < j < L, 

(5.15) QjUh^Uh in n\uji+i. 

Proof. If p e JI{Mj) and p ^ uji+i (open set !), then Kp (t uji+i [Kp e Mj). 
Recall that Kp was deliberately chosen such that there is X G M.i with Kp C K. 
Since Uh is hnear on K, the same holds for Kp and (|5.1|) guarantees 

{QjUh){p) = Uh{p) . 

When restricted to 17 \ Wi+i, the mesh Aij is a refinement of AAh- Hence, agreement 
of the A^?i-piecewise linear function Uh with QjUh in all nodes of M.j outside uii+i 
implies Q,Uh|o\^,+, = 

Consequently, for any Uh G V{JVlh), outside both QiUh and Qi-iUt agree with 

Corollary 5.6. For any G V(A^/i) anrf 1 < ^ < 
supp((Q; - Qi_i)uh) c IJ; . 

In other words, the components of ()5.10p are localized inside refined regions of Q. 
In light of the definition (j4.ip of the refinement zones, we also find 

(5.16) {Qi-Qi_i)uheViMi)\ 

However, having used Q; we cannot expect the splitting to match potential 
homogeneous Dirichlet boundary conditions. This can be remedied using Oswald's 
trick [44, Cor. 30]. We fix Uh G V{Mh) and abbreviate uq = QqUh G V{Ma), 
ui :— {Qi — Qi-i)uh G V{Mi), I > 1. Then, we consider the partial sums 

I 

(5.17) si := ^ ^i-^i) ^ > ■ 

j=o 

Dropping those basis functions in ^v{.Mi) that belong to vertices 'mTo in the rep- 
resentation of s; we arrive at si G V{Mi) G H^^{n). 
Due to Cor. 15. 6[ we observe that 

(5.18) si and s;_i agree on J7 \ uji. 

Hence, away from Ui HT d the same basis contribution are removed from both functions 
when building s; and s/-i, respectively. This permits us to conclude 

(5.19) si and s;_i agree on J7 \ uji. 
Putting it differently, 



(5.20) 



supp(s/ - s/„i) C uji 
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Hence, for all 1 < Z < L we can estimate 

The benefit of zeroing in on lui is that on this subdomain s; has the same "uniform 
scale" hi as ui. Thus, repeated application of uniform L^-stability estimates (|2.9p for 
basis representations and elementary Cauchy-Schwarz inequalities make possible the 
estimates (for arbitrary < e < ^) 

2 

L I 

j=0 j=0 



p6A/'(r,) 



<chi[Y. \\^A\mT,^'' <chi[Y. ^"II"jHl2(.,)'' 
j=i+i j=i+i 

j=i+i 



Here the set M{Yi) comprises the nodes of N{M.i) that lie on a;; n Fd and we make 
heavy use of the geometric decay of hi. The latter also yields 



L L L 

E pi < cE E hf^"^ \\'^o\tmu,) 

1=1 1=1 j=l+l 

= ^EfEv"V7'^"Kii' 



0=2 \l=l 
L 



<Cj2hf\\u,\\l..^^.^<C\u,\l. 



(Q) 



J=2 



by virtue of Lemma 15.41 Except for the last line, all constants depend only on 
and the constants in (|3.3p . Merging the last estimate with (|5.2ip gives us 



(5.22) J2^^'\\'^~'^-^\\lHn)<CWh\li 



(O) 



1=1 



Thus, in light of (|5.20p and the following identity 



So + E('*' " '^'-l) = = SL = Uh, 
1=1 



we have accomplished the proof of the following theorem: 

Theorem 5.7. For any Uh & V{Mh) we can find ui e V{Mi) such that 

L 

(5.23) Uh = E"'' supp(m/) Cuji , 

(=0 
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and 

L 

l"olffi(o) +Zl^r^ ll^illi^(n) <G\uh\\^^^n^ , 
1=1 

with C > independent of L. 

Notice that in combination with the L^-stabiUty (|2.9p of nodal bases and inverse 
inequahties, this theorem asserts an L-uniform estimate of the form (|4.6p for the 
sphtting (|i?TU)l w.r.t. the energy norm H/fi^Qj- From (|5.23p it is clear that the basis 
functions admitted in (j4.10p can represent the functions u; of Thm. 15.71 

Remark 5.8. It is interesting to note that, in contrast to other analyses [1, 11], 
the above proof does not hinge on Assumvtion \3.07J\ Thm. \ 5.7\ remains valid for an 
arbitrary number of hanging nodes on an active edge. Howver, this does not translate 
into asymptotically optimal convergence of local -multigrid in this case, because, in 
order to infer it from Thm. \5. 7[ we also need uniform -stability of the bases. 

5.3. Helmholtz-type decompositions. Helmholtz-type decompositions, also 
called regular decompositions , have emerged as a powerful tool for answering questions 
connected with if (curl, il). In particular, they have paved the way for a rigorous 
multigrid theory for (curl, J7)-elliptic problems [17,25,29,31-33,35,47]. We refer 
to [30, Sect. 2.4] for more information. 

We will need a very general version provided by the following theorem. 

Theorem 5.9. Let meet the requirements stated in Sect. [II Then, for any 
V G i?rr> (curl, ^2), there exists a p £ H^^{Q) and '3> G {H^^{n))^ such that 

(5.24) v = Vp+*, 

< C|| curlv||L2(a), 

where the constant C depends only on fi. 




Fig. 5.2. Buffer zones attached to connected components of (red) Dirichlet boundary part Vn 
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Proof. Given u e /fpo (curl, 51), we define u G iJ(curl, fi), 57 := interior(ri U 
r^i U U . . . ) (see Sect. [1] and Fig. 15.21 for the meaning of fi^), by 



(5.26) u(a;) 



u(a;) for a; G , 

for a; G fli for some i 



Notice that the tangential components of u are continuous across dil, which ensures 
u G £r(curl, Q). Then extend u to u G iJ(curl, R^), see [16]. 

Since curlu G if(div 0, R^), Fourier techniques [24, Sect. 3.3] yield a $ G 
(i?i(M3))3 that fulfills 

(5.27) curl* = curlu, ||*||^i(]g3) < C ||curlu||^2(R3) , 
with C = C{il) > 0. As a consequence 

(5.28) curl(u -*)=0 => u-* = grad q in M.^ . 

On every fli, by definition u = 0, which implies q^Q. G H'^{fli). As the attached 
domains fi^ are well separated Lipschitz domains, see Fig. 15.21 the iJ^-extension of 
q\\j.Q. to g G H^{M.^) is possible. Moreover, it satisfies 



(5.29) ll9llff2(R3) < C'II9|Ih2(u^o,) < C'II*IIhi(r3) < l|curlu||^2 



(n) • 



(5.30) u = — grad q + grad((7 + q) . 
Finally, set 'U' := (* — gradg)|Q, p -.^ q + q, and observe 

(5.31) 11*11 

H^n) - ^ ll'?llff2(E3) < C ||curlu||^2(f2) . 

The constants may depend on 57, Td, and the chosen 57^. □ 

The stable Helmholtz-type decomposition (|5.24p immediately suggests the follow- 
ing idea: when given Vh G \J{Mh), first split it according to ()5.24p and then attack 
both components by the uniformly iJ^-stable local multilevel decompositions explored 
in the previous section. Alas, the idea is flawed, because neither of the terms in (I5.24p 
is guaranteed to be a finite element function, even if this holds for v/j. 

Fortunately, the idea can be mended by building a purely discrete counterpart of 
(|5.24p as in [33, Lemma 5.1] (called there "discrete regular decomposition"). For the 
sake of completeness we elaborate the proof below. 

Lemma 5.10. For any G lJ{Mh), there is G (YiMhlY , Ph G V{Mh), 
and v/i G \J{Aih) such that 

(5.32) Vh=^h + n,,*,, + Wph , 
(5-33) \\Ph\\HHn) <<^\Wh\\H(cur\,n) > 
(5-34) \\^~'^^h\\L2(^n) + \\'^h\\m{n) < C II curl v,,| 1^2(^2) , 

where the constant C depends only on Q, Td, and the shape regularity of Mh- 

Proof, [cf. [33, Lemma 5.1]) We fix a G \J{M.h) and use the stable regular 
decomposition of Thm. 15.91 to split it according to 



(5.35) 



v„ = * + gradp, vi/G (if 1^(57) )3 , p e H^JQ) 
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We have already known that the functions '3/ and p satisfy 

(5.36) < C||curlv,,||^2(n) , llgi"adp||^2(j^) < C ||vh||^(^^ri,n) > 

with constants depending only on and Yd. 

Next, note that in (|5.35p curl* = curlv/j G curlU(A^;i), and, owing to 
Lemma 12.11 II/i'lI' is well defined. Further, a commuting diagram property together 
with Lemma 12.31 implies 

(5.37) curl(/rf-n,0* = ^ 3q ^{Vl) : (/rf - n,0* = grad g . 
The estimate of Lemma 12.11 together with (|5.36l) yields 

(5.38) 

In order to push 'S' into a finite element space, a quasi-interpolation operator 
Qh : {L'^{^)Y ^ V^{^h)Y is the right tool. We simply get it from componentwise 
application of an operator according to Def. I5.1l where any node^element-assignment 
will do. Thus, we can define the terms in the decomposition (|5.32p as 

(5.39) V,, := n,,(* - Q,,*) e U(X,,) , 

(5.40) := Q„* e {V{Mh)f , 

(5.41) gradp/i := grad(p + q) , ph G V{Mh) ■ 

Indeed, grad(p + g) e \J{Mh) such that p + q & V{Mh)- The stability of the decom- 
position (j5.32p can be established as follows: first, make use of Lemma [2. II and (|5.7p 
to obtain, with C = C{pMh) > Oi 

W'^^'^'AIl^q) ^ ll'^"'(^'^"nft)(*-Q/i*)|L2(jj) + ||/i"'(M-Q,0*|L2(n) 

Due to the definition (|5.40p . the next estimate is a simple consequence of (15.61) and 
Thm.EH 

(5.42) W^hWnHn) < CmH^n) < C ||curl v,||^.(^) . 
Finally, the estimates established so far plus the triangle inequality yield 

(5.43) ||gradp„||^.(^) <C||v,|| 

□ 

5.4. Local multilevel splitting of \5{Mh). With the discrete Helmholz-type 
decomposition of Lemma 15.101 at our disposal, we can now tackle its piecewise linear 
and continuous components with Thm. 15.71 

Lemma 5.11. For any G \J{A4h), there exists a constant C only depending on 
the domain, the Dirichlet boundary part Td, the shape regularity of the meshes A4i, 
■Ml, ^ <l < L, and the constants in p.3p . such that 

L 

(5.44) = ^ (v, + Vpi) , V, G Span {53^} , pi e Span {53^-} , 

1=0 
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and 
(5.45) 

L 

l|vo|lL(curl,fl) + IPolHi(0)+Xl^"^ (llvilli^fo) + Ibilli^fo)) <C\\Xh\\ 

Jf(curl,n) ' 

1=1 

where 03^ and 05^ are defined in (|4.9p . 

Proof. We start from the discrete Helmholtz-type decomposition of v^, in (|5.32p : 

V,, = V,, + Uh^h + Vph, *h e (V"(Xh))^ e V{Mh), V,, e U(X,0- 

We apply the result of Thm. [5771 about the existence of stable local multilevel splittings 
of V{A4h) componentwise to >&/,.: this gives 

L 

(5.46) *^ = ^*i, eSpan{Q5i.}' , 

1=0 
L 

(5.47) l*olLi(o)+E'^r'll*'lli^(o) <^l*'^l?^Hn) • 

1=1 

Observe that the functions do not belong to \J{Aii). Thus, we target them with 
edge element interpolation operators 11; onto XJ{Aii), see (|2.7I) . and obtain the split- 
ting described in Lemma [ 



(5.48) ^i=Ui'^i+Vwi , wieV2iMi). 

The gradient terms introduced by (|5.48p are well under control: writing :~ X)^Lo 
the L^-stability of (|5.48p . see Lemma yields 



\\WihHn)<C\\^i\\mn) , 

//i(n) • 

;=o 1=0 1=0 

Because of curlllo'S'o = curl'S'o, we infer from (|5.47p 

L 

(5.49) ||no*oll 

1=1 

Above and throughout the remainder of the proof, constants are independent of L. 

By the projector property 11^011; = Tli, I = 0, . . . , L, and the commuting diagram 
property (|2.1ip . we arrive at 

L 

(5.50) V,, = vh + E^i*; + gradilhSh + ph) , 

1=0 

where 2h is the nodal linear interpolation operator onto V{A4h). Recall (|2.22p to see 
that 

\IhSh+Ph\Hi{n) < C\sh\Hi{n) + \Ph\H^n) < \Wh\\H(curi,n) ■ 
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The local multilevel splitting of X/jS/j + Ph according to Thm. [S771 gives 

L 

(5.51) IhSh+Ph^^Pi, _p/ e Span{<8y} , 

L 

(5.52) bol^i(o) llP'lli^(o) <C\IhSh+Ph\Hii^n} < C' l|v/^llL(curLn) ■ 

1=1 

Still, the contribution V/j does not yet match (j4.1ip . The idea is to distribute V/j 
to the terms 11; 'i'; by scale separation. To that end, we assign a level to each active 
edge of Mh 

(5.53) e{E) mm{£{K) : K e Mh, E is edge oi K} , Ee £{Mh) ■ 
Thus, we distinguish parts of on different levels: given the basis representation 

(5.54) Vft = ashs , {hs} E^^f^^^) = ^ij{Mh) , 

we split 

L 

(5.55) Vft = ^ V; , V; ^ ashs , supp(v;) C uJi . 

1 = EES(M^) 
t{E)=l 

The estimate ||^~^v;i||^2^^j < C ||curlv;i||^2(f^-) from Lemma [5.101 means that v^, is 

"small on fine scales" . Thanks to the L'^-stability (|2.6p of the edge bases, this carries 
over to v;: 

L L 

E'^r'iivdii.(o) <c^Ev' E «liibBiii.(o) 

1=0 1=0 Ee£(Mh),e(E)=l 

L 

(5.56) <CEV' E «l;l|bB|li2(^,) 

1=0 Ee£{Mh),e{E)=i 

L 
1=0 

where G Aiu is coarsest element adjacent to E, cf. (|5.53p . and refinement strips 
are defined by 

(5.57) E; := uji \ aJj^ , <l < L, T^l ■= ^^l , 



see Figs. |4T] and l4?2l 

Yet, in the case of bisection refinement, v; may not be spanned by basis functions 
in Q3u, because the basis functions of JJ{A4h) attached to each edge on S; P| cj;+i, 
< ^ < i do not belong to any *Bu! 

Take any E C f]uJr+T. Let b^;, b^, and b^^ be the basis functions of \J{Mh), 
U(A^;), and U(A^i+i) associated with E, see Fig. 15.31 for a 2D illustration. Denote 
by Ki, . . . , Kn all elements in wj+i and A4i which contain i?, and by . . . , Em their 
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new edges connecting E but not contained in the refinement edges of Ki , . . . , Kn 
(see Fig. 15. 4p . Supposing tfie orientations of each Ei and E point to their common 
endpoint, we have 



(5.58) 



This decomposition is L^-stable with constants merely depending on shape regularity. 




Fig. 5.4. Situation at an edge E lying on the interface between S; and 

Since J2iLi b^^^ G ^' move the component of v; associated with this 

term to v/+i for any E. Then the decomposition (j5.55p and the stability estimate 
(|5.56[) remain valid. 

Summing up, the stability estimate (|5.49p is preserved after replacing with 
Ui^h + 'fi &\JiMi).D 

Eventually, the proof of Thm. 14.21 is readily accomplished. With Lemma 15.111 
at our disposal, we merely appeal to the L^-stabilities expressed in (|2.6p and (|2.9p 
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and inverse inequalities to see that all components in (j5.44p can be split into local 
contributions of basis functions in 03^ and 03^ , respectively. 

6. Quasi-orthogonality. The strengthened Cauchy-Schwartz inequality (|4.7p 
has been established in [61,65] for iJ^-conforming linear Lagrangian finite element 
spaces, in [28, Sect. 6] for H(div)-elliptic variational problems and so-called face ele- 
ments. It is discussed in [31, Sect. 4] for (jl.ip . edge elements, and geometric multigrid 
with global refinement. The considerations for locally refined meshes are fairly similar, 
but will be elaborated for the sake of completeness. 

The trick is, not to consider the one-dimensional spaces spanned by individual 
basis functions as building blocks of the splitting (14. 3p . but larger aggregates. Thus, 
we put the nodal basis functions in *By and *Bu into a small number of classes, such 
that the supports of any two basis functions in the same class do not overlap. Since 
these basis functions are attached to vertices and edges respectively, the definition of 
those classes can be based on a partitioning the vertices/edges of -Mii^j, into disjoint 
sets such that any two vertices/edges of the same set do not belong to the same 
tetrahedron. The is formally stated in the following "colouring lemma" : 

Lemma 6.1. There exist Pj\f, Pg £ N depending only on shape regularity such that 
the sets -f^iMi)^^^ and £{A4i)^^ of vertices and edges of Aii inside the refinement 
zone LJi can he partitioned into subsets 

{p e V{Mi), supp6p c EJ,} = A^^^ U ■ • • U A^,^^ , 
^(■^Olc., := {E e £{Mi), supp&B C cJ,} = 5/ U • • • U , 

and for any K d A4i, K <Z ui, two of its vertices/edges will belong to different subsets. 
Here, bp is the nodal basis function of V{Mi) attached to the vertex p, and b^; is the 
nodal basis function of XJ{A4i) associated with the edge E, see (|2.4p . 

Proof. A crude argument cites the fact that each vertex and each edge belongs to 
only a finite number of elements. A bound for this number can be deduced from the 
shape regularity measure. The rest is elementary combinatorial arguments. □ 

Next, define subspaces of V{A4i) and XJ{A4i) by 

V/ := Span {bp, p G A/"/} C Span {SV} , i = 1, . . . , Pv , 
Uj :=Span{b£;, E e C Span{<8^} , i = 1,...,P£ . 

Note that the basis functions spanning both V^* and are mutually orthogonal 
(w.r.t. a and the if^(ri)-inner product). Thus, it suffices to estabhsh the strength- 
ened Cauchy-Schwarz inequality (|4.7p for the family of subspaces {Hj}^ = {U;}; ■ U 
{grad V[^}i ^ of XJ{Ml)- This will yield the relevant constants in (|4.8p . In other words, 
we analyze the quasi-orthogonality property of the multilevel decomposition 

(6.1) U(A^^) = U(A^o)+EEg^^dl^^ + ^5^Uj . 

1=1 1=1 1=1 i=l 

Note that (|6.ip gives rise to a multigrid algorithm, for which Thm. 14.11 gives exactly 
the same convergence estimate as for the method induced by (|4.1ip ! 

Lemma 6.2. For all e lJ{Mm) and e Uj, 0<m<l<L,i^l,..., Pg, it 
holds that, with C > depending only on the bound for the shape regularity measures 
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of the meshes Mi, 

- — - 

(6.2) (curlv„,curlu;)^2^j-^^ < Chj'h.m'-' ||curl v„i||^2(f^) 1 1 curl U; 1 1 ^^^^ , 

(6-3) (v„,uj)^,(j^^ < Chi W^mWmn) ||curluj||^,(j^^ . 

Proof. Pick any (open) tetrahedron K e Mm- Use the basis representation of U; 
to isolate "interior" and "boundary" parts 



where 



U; 



l,bd 



J2 ( / ui • ds) • bB and uj^^^^ = ^ ( / uj • ds) • b^; . 

EcdK.EeSf ^ EcK.EeSf ^ 



Since curl is a constant vector in K and uj x n = on dK, by Green's formula, 
it is easy to see 

/ curl U; • curl Vm da; = / curl U; ;,^ • curl da; = / curl U;^^ ■ curl v^dx, 

Jk Jk ' is 

where 

E := |J{suppbB : EcdK,E€ 

is contained in a narrow strip along the boundary of K of width « hi. Hence, we 
arrive at the area ratio 

\n<chih-^\K\ . 

Here and throughout the remainder of the proof, C > depends on shape regularity 
only. Thus, using the Cauchy-Schwartz inequality and noting that the basis functions 
in U; are mutually orthogonal, we have 

(6.4) j curlu; -curlvrnda; < ||curluj b^H^^^^^ [Sj^/^ Icurlv^l 

||curlui||^,(^^ \K\^/^ \cuvWm\ 

= C\/t^ ||curluj||^2(^) ||curlv„||^2(x) . 

To estimate the i^-inncr product, we recall the following simple fact about the norms 
of edge basis functions on level I: 

||b||^.(^)<C/in|curlb||^.(^) yhe^{j. 
Since the basis functions of U; do not interact, we have 

(6-5) j^ul ■ Vmdx < ||uj||^,(^j \\-Vm\\L2(K) < ^hi || CUrl uj || ^.^^^ l|Vm||z,2(^) • 
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Now (|6.2p and (|6.3p follows by summation over all elements of Aim and another 
Cauchy-Schwarz inequality. □ 

After replacing with in the proof of Lemma [6.2[ similar arguments establish 
the following estimate: 

Lemma 6.3. For all e U(A^,„) and u'l g 0<m<l<L,i=l,..., Pj^, it 
holds that, with C > depending only on the hound for the shape regularity measures 
of the meshes A4i, 

(6.6) (v„,gradu^)^2(^^ < Chfhm" ||v,„||^2(o) ||gi"ad uj||^2(n) ■ 



Proof. Again, pick K G A^m- By separating interior and boundary parts of u] as 
above and noting divv„i|^ = on K, we find by Green's formula 



v,„ ■ grad u^dx — / • grad dx . 
K Js 



As above, we infer 



K 




v,„ • grad -uj da; < C — ||v„||^2(^) ||gradMj||^2(^) 



Summation over all K and a Cauchy-Schwarz inequality finish the proof. □ 

Because of the geometric decay of the meshwidths hi of the (uniformly refined) 

meshes A4i, these estimates clearly imply the desired quasi-orthogonality for (|6.ip . 
Theorem 6.4. (Strengthened Cauchy- Schwartz inequality) For any u] £ UJ or 

u| G gradVJ' and any vj G or vj G gradVJ-', < i,j < Ps or < i, j < Pj^ , 

resp., the estimate 

(6.7) a(ui,v^„)<C0l'-"l/2||uj||^||v^„||^ < j < A^,, < Z, m < i 

holds, where C — C(maxi p;) > and Q < 6 < 1 is the decrease rate of the meshwidths 
defined in 



7. Numerical experiments. In the reported numerical experiments the im- 
plementation of adaptive mesh refinement was based on the adaptive finite element 
package ALBERTA [48], which uses the bisection strategy of [36], see Sect. [31 

Let A^o be an initial mesh satisfying the two assumptions (Al) and (A2) in [36, 
P. 282], the adaptive mesh refinements are governed by a residual based a posteriori 
error estimator. In the experiments we assume the current density f G (div, $7) and 
use the estimator given by [17, §5]: given a finite element approximation \ih G U(A^/i), 
for any T e Mh 



FcdT 



where F is a face of T, v is the unit normal of F, and [uhlr is the jump of u;i across 
F. The global a posteriori error estimate and the maximal estimated element error 
on Aifi are defined by 




(7.1) Vh-= \ Vt] > r?max = ^max^ ?7T- 
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Using 77/j and rfoiaxj we use [17, Algorithm 5.1] to mark and refine Adh adaptively. 

In the following, we report two numerical experiments to demonstrate the compet- 
itive behavior of the local multigrid method and to validate our convergence theory. 

EXAMPLE 7.1. We consider the Maxwell equation on the three-dimensional "L- 
shaped" domain V, = (-1, 1)^ \ {(0, 1) x (-1, 0) x (-1, 1)}. The Dirichlet boundary 
condition and the righthand side f are chosen so that the exact solution is 

u := v|ri/2 sin(<?!)/2)| 

in cylindrical coordinates {r,(j),z). 

Table 17.11 shows the numbers of multigrid iterations required to reduce the initial 
residual by a factor 10"* on different levels. We observe that the multigrid algorithm 
converges in almost the same small number of steps, though the number of elements 
varies from 156 to 100,420. 

Table 7.1 

The number of adaptive iterations N^, the number of elements N^i, the number of multigrid 
iterations Nurs required to reduce the initial residual by a factor 10~®, the relative error between the 
true solution u and the discrete solution u^; iJ^cl = 11^ ~ '^hllH(curl n) / ll"llH{curl n) (Example 

[HP- 



iVit 


2 


5 


10 


15 


20 


25 


30 


35 


iVel 


156 


388 


1,900 


4,356 


9608 


19,424 


48,088 


100,420 


Erel 


0.4510 


0.3437 


0.2456 


0.1919 


0.1600 


0.1350 


0.1094 


0.0915 




11 


21 


19 


19 


19 


19 


19 


19 



Fig. 17.11 (left) plots the CPU time versus the number of degrees of freedom on 
different adaptive meshes. It shows that the CPU time of solving the algebraic system 
increases roughly linearly with respect to the number of elements. Fig. 17.11 (rught) de- 
picts a locally refined mesh of 100,420 elements created by the adaptive finite element 
algorithm. 




Fig. 7.1. Example ]?. 1\ left: execution time for local multigrid method, right: instance of a locally 
refined mesh (100,420 elements) 

EXAMPLE 7.2. This example uses the same solution as Example \ 7. 1\ 

u v|r^/2 sin(0/2)| 
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in cylindrical coordinates (r, 0, z). But the computational domain is changed to a three- 
dimensional non-Lipschitz domain with an inner crack-type boundary, which is defined 
by 

n = (-1, 1)3 \ {(x, 0, z) : < X < 1, -1 < z < 1}. 

The Dirichlet boundary condition and the source function f are the same as above. 

Table [721 records the numbers of multigrid iterations required to reduce the initial 
residual by a factor 10~* on different levels. We observe that the multigrid algorithm 
converges in less than 30 steps, with the number of elements soaring from 128 to 
135,876. 

Table 7.2 

The number of adaptive iterations Nit, ihe number of elements N^i, the number of multigrid 
iterations Nurs required to reduce the initial residual by a factor 10~*, the relative error between the 
true solution u and the discrete solution u^.- i?rcl = l|u ~ '^'j|lH(curl fi) / ll'^llH(curl fl) (Example 



Nit 


2 


5 


10 


15 


20 


25 


30 


33 




128 


404 


1,236 


3,416 


12,420 


29,428 


81,508 


135,876 


Ercl 


0.4616 


0.3762 


0.2992 


0.2347 


0.1752 


0.1394 


0.1095 


0.0958 


Nitrs 


14 


30 


25 


26 


26 


27 


27 


27 



Fig. 17.21 (left) shows the CPU time versus the number of degrees of freedom on 
different adaptive meshes. Obviously, the CPU time for solving the algebraic system 
increases nearly linearly with respect to the number of elements. 

Fig. 17.21 fright') displays a locally refined mesh of 135,876 elements using adaptive 
finite element algorithm. In addition, the restriction of the mesh to the cross-section 
{y = 0}, which contains the inner boundary, is drawn. This reveals strong local 
refinement. 

This experiment bears out that the local multigrid is also efficient for the problems 
in non-Lipschitz doamins, which are outside the scope of our theory. 




Fig. 7.2. Example \ 7.2\ left: CPU time for solving the algebraic system by multigrid method, 
right: a locally refined mesh (135,876 elements) 
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